DESY-09-122 
SFB/CPP-09-70 



S quark pair production at the LHC 

U. Langenfeld 

Deutsches Elektronensynchrotron DESY, Platanenallee 6, D-15738 Zeuthen 

Abstract. We present NNLO cross sections for squark- antisquark production at the LHC. We have calculated new analytic 
expressions for the scale dependent scaling functions at one and two loop. 

Keywords: Collider physics, supersymmetric particles, pair production, higher order calculation 
PACS: 12.38.-t, 12.38.Bx, 12.60.Jv, 14.80.Ly 



1. SQUARK PAIR PRODUCTION CROSS SECTION AT THE LHC 

If Supersymmetry is realised in Nature then it is expected that squark and gluino pairs are produced in large numbers 
at the LHC. It is possible to probe masses up to the TeV range. Squarks are assumed to be heavier than w 400 GeV IH 
so these particles are produced near the kinematical production threshold. Therefore one can use the same methods to 
calculate higher order cross sections as developed for tt production 12j,|3[]. The partonic LO and NLO cross sections 
are known for long times 0]. Approximate NNLO corrections have been calculated in OJ]. The LO partonic cross 
section and the NLO theshold expansion are known analytically U In this article we present analytical formulae 
for the scale dependence determining NLO scaling functions and the threshold expansion of the scale dependence 
determining NNLO scaling functions. For related work, see also Ref. JfJ. In Ref. 0], the soft anomalous dimension 
has been calculated to NNLO accuracy. 

The partonic cross section s with identified renormalisation and factorisation scale can be expanded as 
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L M = log(ii 2 /m 2 q ), (1) 



where ij denote the initial states gluon - gluon or quark - antiquark. The full dependence on the renormalisation and 
factorisation scale is the same as for ff-production, see Ref. |81|. The hadronic cross section is given as a convolution 
of the partonic cross section with the parton luminosities L,-,-: 
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<y pp ^qq*x(s,mq,mg)= ^ / dsLij &ij(s,mq,mg,ii). 
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We performed a scan of the LO, NLO, and NNLO squark pair production cross section in the nig - trig - plane 
using Prospino Jit] and the formulae presented in Q], see Fig.QJa) -(c). One clearly sees the strong enhancement 
of the NLO and NNLO cross section compared to the LO cross section: For squarks and a gluino with mass 200 GeV 
and 250 GeV, respectively, we have a LO cross section of about 500 pb, but about lOOOpb at NLO and NNLO. At 
NNLO, the lOOOpb region is even enlarged to higher gluino masses. For a squark mass of 400GeV and a gluino 
mass of 500 GeV, we find for the LO, NLO, and NNLO cross section 28.9pb, 43.1 pb, and 46.9pb, respectively. The 
NNLO cross section is 9% larger than the NLO cross section. The contour lines of constant cross section are running 
nearly parallel to the gluino mass axis: The cross section shows a rather mild dependence on the gluino mass. There 
is a weak enhancement of the cross section for m q = nig as one can see from the small bump at nig = mq. The cross 
section decreases for more than three orders of magnitude for squark masses from 200 — 1000 GeV. This strong mass 
dependence is well-known from hadronic tt pair production. 

In Fig. [TJd) we show the full jij - jj, r scale dependence of the NNLO cross section for the example point mq = 

400 GeV, nig = 500GeV. The scale uncertainty is about —8% for (jU.y,jU r ) = (^mq,2m^) and about +8% for (jU.y,/J r ) = 
(2niq, \niq). This is considerably larger than the usual scale uncertainty taken at jx r = jj.^ = jX (in our example w —4% 
at jU = 1/2 and + 1 % at ju = 2). This shows that a full treatment of the scale dependence leads to more reliable estimates 
of the scale uncertainty. 





2. ANALYTICAL FORMULAE 

In this section, we present analytical formulae for the /y scaling functions and the threshold logarithms of the 

(22) (11) 

and f>. scaling functions. The f\ scaling functions are determined by the renormalisation group equation: 

f (11) - 1 (R f (00) P {Q) 6?,f m \ ii-oo nn (Vi 

The are the leading order splitting functions, see Ref. IToll . ® denotes the standard Mellin convolution. The scale 
dependent NLO scaling function Eq. (0) depends only on LO functions. Performing the integrations yields as new 
analytic results the Eq. (0]i and (O. 
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L 4 = Li 2 (^)-Li 2 (±±£), 
L 6 = log 2 (1+/3)- log 2 (1-/3). 
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The leading order scaling functions /gg '' and /i2 can be found in Ref. The NNLO scale dependent scaling 
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functions follow from the RGE relations 
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i,j,k,l,n are parton indices with implied summation over repeated indices. The threshold expansion is derived by 
computing the Mellin transformation of each of the involved factors and inverting the products back to p space. The 
scaling function fgq is very small near threshold so we did not include them. The constants af ,us can be found in 



Ref. |a. The coefficients of the QCD J3 - function are given as J3 = 11 - (2/3 )«/ and J3i = 102 - (38/3)«/. We used 
the threshold expansion to fit the numerical determined values of the NNLO scaling functions. 
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